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Abstract 



Lie symmetry group method is applied to study Newtonian incompressible fluid's equations flow in turbulent boundary 
^ ' , layers. The symmetry group and its optimal system are given, and group invariant solutions associated to the symmetries are 



obtained. Finally the structure of the Lie algebra symmetries is determined. 
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. 1 Introduction 
> 

l/~j In physics and fluid mechanics, a boundary layer is that layer of fluid in the immediate vicinity of a bounding surface. 

^^ ■ In the Earth's atmosphere, the planetary boundary layer is the air layer near the ground affected by diurnal heat, 

^^ ' moisture or momentum transfer to or from the surface. On an aircraft wing the boundary layer is the part of the flow 

r — I , close to the wing. The boundary layer effect occurs at the field region in which all changes occur in the flow pattern. 

f^ ' The boundary layer distorts surrounding non- viscous flow. It is a phenomenon of viscous forces. This effect is related 

f^ I to the Reynolds number (In fluid mechanics and heat transfer, the Reynold's number is a dimensionless number 

that gives a measure of the ratio of incrtial forces to viscous and, consequently, it quantifies the relative importance 

of these two types of forces for given flow conditions). Laminar boundary layers come in various forms and can be 

loosely classified according to their structure and the circumstances under which they are created. The thin shear 

S^ I layer which develops on an oscillating body is an example of a Stokes boundary layer, whilst the Blasius boundary 

\^ ' layer refers to the well-known similarity solution for the steady boundary layer attached to a flat plate held in an 

5t i oncoming unidirectional flow. When a fluid rotates, viscous forces may be balanced by the Coriolis effect, rather 

than convective inertia, leading to the formation of an Ekman layer. Thermal boundary layers also exist in heat 

transfer. Multiple types of boundary layers can coexist near a surface simultaneously. The deduction of the boundary 

layer equations was perhaps one of the most important advances in fluid dynamics. Using an order of magnitude 

analysis, the well-known governing Navier— Stokes equations of viscous fluid flow can be greatly simplified within 

the boundary layer. Notably, the characteristic of the partial differential equations (PDE) becomes parabolic, rather 

than the elliptical form of the full Navier— Stokes equations. This greatly simplifies the solution of the equations. By 

making the boundary layer approximation, the flow is divided into an inviscid portion (which is easy to solve by a 

number of methods) and the boundary layer, which is governed by an easier to solve PDE. 

Flow and heat transfer of an incompressible viscous fluid over a stretching sheet appear in several manufacturing 
processes of industry such as the extrusion of polymers, the cooling of metallic plates, the aerodynamic extrusion 
of plastic sheets, etc. In the glass industry, blowing, floating or spinning of fibres are processes, which involve the 
flow due to a stretching surface. Mahapatra and Gupta studied the steady two-dimensional stagnation-point flow 
of an incompressible viscous fluid over a flat deformable sheet when the sheet is stretched in its own plane with a 
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velocity proportional to the distance from the stagnation-point. They concluded that, for a fluid of small kinematic 
viscosity, a boundary layer is formed when the stretching velocity is less than the free stream velocity and an inverted 
boundary layer is formed when the stretching velocity exceeds the free stream velocity. Temperature distribution in 
the boundary layer is determined when the surface is held at constant temperature giving the so called surface heat 
flux. In their analysis, they used the finite-differences scheme along with the Thomas algorithm to solve the resulting 
system of ordinary differential equations. 

The treatment of turbulent boundary layers is far more difficult due to the time-dependent variation of the 
flow properties. One of the most widely used techniques in which turbulent flows are tackled is to apply Reynolds 
decomposition. Here the instantaneous flow properties are decomposed into a mean and fluctuating component. 
Applying this technique to the boundary layer equations gives the full turbulent boundary layer equations 



(1) 



where p is the density, p is the pressure, v is the kinematic viscosity of the fluid at a point and u and v are average of 
the velocity components in Reynold decomposition. Here u' and v' are the velocity fluctations such that; u = u + u' 
and v = V + v' . By using the scale analysis (a powerful tool used in the mathematical sciences for the simpliflcation 
of equations with many terms), it can be shown that the system (1) reduce to the classical form 
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The term u'v' in the system (2) called Reynolds shear stress, a tensor that conventionally written 



Rij=pu\u'y (3) 

The divergence of this stress is the force density on the fluid due to the turbulent fluctuations. Using NavierStokes 
equations for a fluid whose stress versus rate of strain curve is linear and passes through the origin (Newtonian fluid) 
the tensor (3) reduces to 



Rij = p 



duj 
dxi 



^3 

where p is the fluid viscosity. 

This paper is concerned with the symmetry of Newtonian incompressible fluid's equations flow in turbulent 
boundary layers of the form (2). Lie-group theory is applied to the equations of motion for determining symmetry 
reductions of partial differential equations. The solution of the turbulent boundary layer equations therefore necessi- 
tates the use of a turbulence model, which aims to express the Reynolds shear stress in terms of known flow variables 
or derivatives. The lack of accuracy and generality of such models is a major obstacle in the successful prediction of 
turbulent flow properties in modern fluid dynamics 

2 Lie Symmetries of the Equations 

A PDE with p— independent and g— dependent variables has a Lie point transformations 

Xi — Xi + e^i{x, u) + O(e^), Ua — Ua + eipa{x, u) + O(e^), 
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where C* = -rr 
oe 



for i — 1 , . . . , p and ipa — —;:: — 
e=o oe 



by its associated infinitesimal generator 



for a = 1, ...,q. The action of the Lie group can be considered 
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^l{x,u)--— + ^LPa{x,u)- (4) 
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on the total space of PDE (the space containing independent and dependent variables). Furthermore, the character- 
istic of the vector field (4) is given by 

P C) a 

i—1 

and its n— th prolongation is determined by 

4=1 * a=l ttJ=i=0 -^ 

where (f^ — DjQ"" + X^f^i ^i''^ji- i^J i^ the total derivative operator describes in (5)). 

Let us consider a one-parameter Lie group of infinitesimal transformations (x, y, u, v,p) given by 

x = x + e^i{x,y,u,v,p) + 0{e^), y = y + e£,2(x,y,u,v,p) + 0{e'^), 

u = u + erii{x,y,u,v,p) + C'(e^), v ^ v + eri2{x,y,u,v,p) + 0{e'^), 

p^p + ei]3{x,y,u,v,p) + 0{e^), 

where e is the group parameter. Then one requires that this transformations leaves invariant the set of solutions 
of the system (2). This yields to the linear system of equations for the infinitesimals ^i{x,y,u,v,p),^2{x,y,u,v,p), 
rii{x,y,u,v,p),ri2{x,y,u,v,p) and ri3{x,y,u,v,p). The Lie algebra of infinitesimal symmetries is the set of vector 
fields in the form of 

. d ^ d d d d 
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This vector field has the second prolongation 
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with the coefficients 

Vl = Dx{ipi - ilUx - £,2Uy) + ilUxx + ^2Uxy, V>\ = Dy(ipi - S,iUx ~ S,2Uy) + (lUxy + ^2Uyy, 

VT = Dl{ipi - S,lUx - C2Uy) + S.lUxxx + C2Uxxy, ^V ^ DxDy{(pi - ^lUx - S,2Uy) + Cl^xxy + S,2Uxxy, 

ipf == D^lipi - ^lUx - S,2Uy) + S,lUxxy + S,2Uxyy, fl = Dx{(p2 - ^l^x " 6^y) + ^l^xx + S.2Vxy, 

^2 = Dy{ip2 - ^iVx - i,2Vy) + ^iVxy + S.2Vyy, (^2^^ = ^x ('/'2 - ^l^x " ^2Vy) + £,lVxxx + i2Vxxy, 

ifil^ = DxDy{ip2 - ^iVx ~ i,2Vy) + S.lVxxy + S,2Vxxy, ^^2 = ^%{V2 - S.lVx - ^2Vy) + ^iVxxy + ^2Vxyy, 

ifl = Dx{^?, - ilPx - S.2Py) + S.lPxx + S,2Pxy, (fil = Dy{(fs - S^ipx - £,2Py) + S,lPxy + S,2Pyy, 

VT = Dl{ip3 - S^iPx - S,2Py) + S.lPxxx + S,2Pxxy, fl^ = DxDy{(p3 - ^ipx - £,2py) + £,lPxxy + £,2Pxxy, 
fV = DI(lP3 - ^ipx - ^2py) + £.lPxxy + £.2pxyy, 



where the operators D^ and Dy denote the total derivative with respect to x and y: 

d _ d _ d _ d ^ d _ d _ d _ d 

Dx = TT- +UxTrr +Vx—: +PxTrz ••• , Dy = —- + Uy—r + Vy—z + Py^r:! ■■■ ■ 
ax ou ov op oy ou ov op 



(5) 



Using the invariance condition, i.e., applying the second prolongation v'^', to system (2), and by solving the linear 
system 
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the following system of 17 determining equations yields: 

ily = 0, 6fi = 0, Cltj = 0: Cls = 0: 

^^xx = 0; 62; = 0, ^2fi == 0, ^2s = 0, 

6p = 0, £,2yy = 0, 771^ = 0, -qiy = 0, 

m = -S.2y 

The solution of the above system gives the following coefficients of the vector field v: 

Ci = Ci + C4X, 6 = C'2 + C^y, m = Ciu - 2C5U, m = -Cgu, % = C'a + 2Cip - 4C5P, 

where C'i for i ~ 1,...,5 are arbitrary constants, thus the Lie algebra g of the Newtonian incompressible fiuid's 
equations flow in turbulent boundary layers is spanned by the five vector fields 
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V5 = y-R 2u-— - V-— - 4p— , 

dy du dv dp 

The commutation relations between these vector fields is given by the table (1), where entry in row i and column j 
representing [v^ , v^] . 

The one-parameter groups Gi generated by the base of g are given in the following table. 
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Since each group Gi is a symmetry group and if u = f{x, y),v = g{x, y) and p = h{x, y) are solutions of the system 
(2), so are the functions 

ui^ f{x + e,y), ui= f{x,y + e), ui^f{x,y), ui ^ e'" /{xe" ,y), u = e"^" f{x,ye''), 

vi= g{x + e,y), vi= g{x,y + e), vi^g{x,y), vi^ g{xe^,y), v ^ e^g{x,ye^), 

pi =h{x + e,y), pi =h{x,y + e), pi = h{x,y) +e, pi ^ e^^^hixe" ,y), p = e^^h{x,ye^), 

where e is a real number. Here we can find the general group of the symmetries by considering a general linear 
combination CiVi + • • • + C1V5 of the given vector fields. In particular if g is the action of the symmetry group near 
the identity, it can be represented in the form g — exp(£V5) o • • • o exp(£Vi). Consequently, if u = f{x, y), v — g{x, y) 
and p = h(x, y) be a solution for system (2), so is 

u^e^fUx + e)e^, {y + e)e^], 

w = g((a; + £)e",(y + £)e"), (6) 

p^e'^''h({x + e)e'',{y + e)eA +£6^", 

a fact that can, of course, be checked directly. 



3 Symmetry reduction for Newtonian incompressible fluid's equations flow in turbulent boundary 
layers 

Lie-group method is applicable to both linear and non-linear partial differential equations, which leads to similarity 
variables that may be used to reduce the number of independent variables in partial differential equations. By 
determining the transformation group under which a given partial differential equation is invariant, we can obtain 
information about the invariants and symmetries of that equation. 

The first advantage of symmetry group method is to construct new solutions from known solutions. Neither 
the first advantage nor the second will be investigated here, but symmetry group method will be applied to the 
Eq. (2) to be connected directly to some order differential equations. To do this, a particular linear combinations 
of infinitesimals are considered and their corresponding invariants are determined. The Newtonian incompressible 
fiuid's equations flow in turbulent boundary layers expressed in the coordinates {x,y), so to reduce this equation 
is to search for its form in specific coordinates. Those coordinates will be constructed by searching for independent 
invariants (r, s) corresponding to an infinitesimal generator. So using the chain rule, the expression of the equation 
in the new coordinate allows us to the reduced equation. Here we will obtain some invariant solutions with respect to 
symmetries in the way of similarity transformations to leave invariant the system (2) under the group of symmetries. 
First we obtain the similarity variables for each term of the Lie algebra g, then we use this method to reduced the 
PDE and find the invariant solutions. Here our system has two-independent and four-dependent variables, thus, the 
similarity transformations and invariant transformations including of invariants up to order 2 comining in table (2) 
and (3), where * and ** are refer to ordinary invariants and first order differential invariants respectively. 



Table 2 

Reduction of system 



vector field 



symmetry transformations 



similarities 



vi r = X + e,s = y, f{r, s) — u, g{r, s) = v, h{r, s) — p 

V2 r — x,s — y + e, f{r, s) — u, g{r, s) — v, h{r, s) = p 

V3 r — x,s — y, f{r, s) — u, g{r, s) — v, h{r, s) = p + e 

V4 r — xe^ , s = y, f{r, s) = e^'^u, g{r, s) — v, h{r, s) 
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e-2-p 
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X — e'^ ,y — r,u — e" f{r),v = g{r),p — e^'^h{r) 



r — x,s — ye^ , f{r, s) = e ^u, g(r, s) = e^v, h{r, s) = e ^p 
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Table 3 
Invariants 



vector field ordinary invariant 



1st order diff. invariant 



2nd order diff. invariant 



Vl 
V2 
V3 
V4 

V6 



y,u,v,p 
x,u, v,p 
x,y,u,v 
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*, Ux, XVx, —, —, Vy, -^ 
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**, Uxx^ Vxx^Pxx^ Uxy^ Vxy^Pxy: '^yy ^ '^yy ^ Pyy 
**, Uxx^ Vxx^Pxx^ Uxy^ Vxy^Pxy: '^yy i '^yy ^ Pyy 
**, Uxx^ Vxx^Pxx^ Uxy^ Vxy^Pxy: ^yy^'^yy^Pyy 

y,^,v,p,Ux,xVx,^,^,Vy,^, 

2 Pxy '^ yy "Pyy 

XUxx 5 X '^xx 5 Pxx 5 Uxy •) XVxy -, ^ •) ~ 5 "^yy ; „2 

*, y'^ux,yvx,y^px,y^uy, y'^Vy, y^Vy, 
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y ^yyj y Pyy 



4 Characterization of invariant system of differential equations 

Differential invariants help us to find general systems of differential equations virhich admit a prescribed symmetry 
group. One say, if G is a symmetry group for a system of PDEs with functionally differential invariants, then, the 
system can be rewritten in terms of differential invariants. For finding the differential invariants of the system (2) 
up to order 2, we should solve the following systems of PDEs: 
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where / is a smooth function of (x, y, u, v,p), 
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The solutions of (7), (8) and (9) are 



/ = const. 



J —(hi— a "^^ P^ ^y '"y Py \ /-|n\ 
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thus, the most general first and second order system (2) admitting (6) as a symmetry group can be rewritten with 
two arbitrary functions (10) and (11). Characterization of invariant system of differential equations is an alternative 
method for reduction of differential equations invariant under a symmetry group. 



5 Optimal system of steady two-dimensional boundary-layer stagnation-point flow equations 

Let a system of differential equation A admitting the symmetry Lie group G,be given. Now G operates on the set of 
solutions S of A. Let s-G be the orbit of s, and iJ be a subgroup of G. Invariant iJ— solutions s G S are characterized 
by equality s ■ S = {s}. If ft, G G is a transformation and s G 5, then h- {s ■ H) = (h- s) ■ (hHh^^). Consequently,every 
invariant _ff— solution s transforms into an invariant /i_ff/i^^— solution (Proposition 3.6 of [5]). 

Therefore, different invariant solutions are found from similar subgroups of G. Thus, classification of invariant 
_ff —solutions is reduced to the problem of classification of subgroups of G,up to similarity. An optimal system of 
s— dimensional subgroups of G is a list of conjugacy inequivalent s— dimensional subgroups of G with the property 
that any other subgroup is conjugate to precisely one subgroup in the list. Similarly, a list of s— dimensional subal- 
gebras forms an optimal system if every s— dimensional subalgebra of Q is equivalent to a unique member of the list 
under some element of the adjoint representation: 1) — Ad(g) • f). 

Let H and H be connected, s— dimensional Lie subgroups of the Lie group G with corresponding Lie subalgebras 
() and 1) of the Lie algebra g of G. Then H — gHg~^ are conjugate subgroups if and only t) — Ad{g) ■ 1) are conjugate 
subalgebras (Proposition 3.7 of [5]). Thus, the problem of finding an optimal system of subgroups is equivalent to 
that of finding an optimal system of subalgebras, and so we concentrate on it. 



5.1 One- dimensional optimal system 

For one-dimensional subalgebras, the classification problem is essentially the same as the problem of classifying the 
orbits of the adjoint representation, since each one-dimensional subalgebra is determined by a nonzero vector in Lie 
algebra symmetries of steady two-dimensional boundary-layer stagnation-point flow equations and so to "simplify" 
it as much as possible. 

The adjoint action is given by the Lie series 

Ad(exp(ev,)vj) = Vj - e[vj, v^-] + — [v^, [v^, v^-]] , (12) 

where [vi,Vj] is the commutator for the Lie algebra, e is a parameter, and i,j = 1, • • • ,5. Let -Pf : g — > g defined 
by V !->• Ad(cxp(evi)v) is a linear map, for i = 1, • • • , 5. The matrices M? of F?, i = 1, • • • ,5, with respect to basis 



{vi,--- ,V5} are 
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/ 1 o\ 
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-s 1 

V 1/ 

/ e= \ 
10 
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10 

\ 1/ 



Ml 



Mi 



( 1 o\ 
10 
10 
10 

V -£ 1 / 

/l o\ 
e= 
e-*^ 
10 

VO 1/ 



MS 



/l o\ 
10 
10 
-2e 1 

V 4e 1/ 



(13) 



by acting these matrices on a vector field v alternatively we can show that a one-dimensional optimal system of q is 
given by 



aiV2 + a2V3, X4 — aivi + a2^2 + asvs, 



(14) 



Xi=V3, X2 = aiVi +a2V2, X3 

where a^'s are real constants. 

5.^ Two-dimensional optimal system, 

Next step is to construct two-dimensional optimal system, i.e., classification of two-dimensional subalgebras of q. 
The process is by selecting one of the vector fields in (14), say, any vector field of (14). Let us consider Xi (or 
Xi,i = 2,3,4,5). Corresponding to it, a vector field X — aiVi + • ■ • + 05 V5, where a^'s are smooth functions of 
{x,y,u,v,p) is chosen, so we must have 



[Xi,X] ^XXi+fiX, 
the equation (15) leads us to the system 

Cji^ajUk = \ai + fiai (i = 1, 2, 3, 4, 5). 



(15) 



(16) 



The solutions of the system (16), give one of the two-dimensional generator and the second generator is Xi or, 
Xi, i — 2,3,4,5 if selected. After the construction of all two-dimensional subalgebras, for every vector fields of (14), 
they need to be simplified by the action of (13) in the manner analogous to the way of one-dimensional optimal 
system. 

Consequently the two-dimensional optimal system of g has three classes of g's members combinations such as 



Vl,V2, 
Vl,V4, 
V3,V5, 



Vl,V3, V2,V3, 
Vl,V5, V2,V4 



V3,V4, 
V2,V5, 



V4,V5, Vi,Yfi^2 1^^^^' /^iVi +^2V2,V3 -h/33(v4-^V5), 



(17) 



/3lV2 + /32V3, Vi + |/33(V4 + V5). 



5.3 Three-dimensional optim,al system. 

This system can be developed by the method of expansion of two-dimensional optimal system. For this take any 
two-dimensional subalgebras of (17), let us consider the first two vector fields of (17), and call them Yi and Y2, thus, 
we have a subalgebra with basis {Yi, 1^2}, find a vector field Y — aivi -I- • • • -|- asVs, where Ui's are smooth functions 



of {x,y,u,v,p), such the triple {Yi,F2,y} generates a basis of a three-dimensional algebra. For that it is necessary 
an sufficient that the vector field Y satisfies the equations 

[Yi,Y]=XiY + ^iiYi + i^iY2, 

[Y2,Y]=\2Y + ^l2Yl+U2Y2, (18) 

and following from (18), we obtain the system 

C]k|iiak = \la, + ^llPl + VlPl, a = 1,2,3,4,5 (19) 

Cife/3fafe = A2ai + /i2,3; + t^2^s, a^ 1,2,3,4,5. 

The solutions of system (19) is linearly independent of {Yi^ 1^2} and give a three-dimensional subalgebra. This process 
is used for the another two couple vector fields of (17). 

Consequently the three-dimensional optimal system of q is given by 

Vl,V2,V3, Vl,V2,V4, Vl,V2,V5, Vl,V3,V4, Vl,V3,V5, 

(zUj 

Vi,V4,V5, V2,V3,V4, V2,V3,V5, V2,V4,V5, V3,V4,V5. 



5.^ Four- dimensional optimal system 

Four-dimensional optimal system can be developed by the method of expansion of three-dimensional optimal system 
analogous to three-dimensional. First step is selecting a three-dimensional subalgebra of (20), and call them Zx^Z^ 
and Z3, then find a vector field Z — a\W\ + • ■ • + 05 V5, where a,'s are smooth functions of {x,y,u,v,p), such the 
foursome {Zi, Z2, Z^, Z} generates a basis of a four-dimensional algebra. It is necessary an sufficient that the vector 
field Z satisfies the equations 

1Z2 +71Z3, 

(21) 



ZuZ]-- 


= XiZ + fiiZi + vxZ2 + 7i2'3 


Zi,Z] = 


= \2Z + [i2Z\ + 1^2 ^2 + 72^3 


Zi,Z] = 


= \zZ + /isZi + 1/3^2 + 73^3 



following from (21), we obtain the system 



fcafc = AiCi +,(ii H- t^i -1-71, a = 1,2,3,4, 5, 

feOfc^ A2ai +,U2 + 1^2 + 72, a = 1,2,3,4, 5, (22) 

feCfe^Asai + /i3 + 1^3 + 73, a =1,2,3,4,5. 



The solutions of system (22) is linearly independent of {Zi, Z2,Z3} and give a four-dimensional subalgebra. This 
process is used for the another triple vector fields of (20). All previous calculations lead to the table (4) for the 
optimal system of g. 

6 Lie Algebra Structure 

Q has a no any non-trivial Levi decomposition in the form of g = r K gi, because q has no any non-trivial radical, 
i.e., if r be the radical of g, then g = r. 

If we want to integration an involuting distribution, the process decomposes into two steps: 

• integration of the evolutive distribution with symmetry Lie algebra g/r, and 

• integration on integral manifolds with symmetry algebra r. 



Table 4 

Optimal system of subalgebras 



dimension 



1 



(V3) 

subalgebras (aivi + 02 V2) 
(aiV2 + a2-V:i) 
(aivi + a2V2 + asva) 



(V1,V2) (vi,V2,V3 

(V1,V3) (V1,V2,V4 

(V2,V3) (vi,V2,V5 

(V2,V4) (vi,V3,V4 

(V3,V4) (V1,V3,V5 

(V1,V4) (vi,V4,V5 

(vijVs) (V2,V3,V4 

(V2,V5) (V2,V3,V5 

(V3,V5) (V2,V4,V5 

(V4,V5} (V3,V4,V5 
(vi,EL2Av,> 

(/?lVl +/?2V2,V3 +/?3(V4 + V5)) 
(/?1V2 + /?2V3, VI + |;53(V4 + V5)) 



(vi,V2, V3, V4 
(vi,V2, V3,V5 
(vi,V2, V4,V5 
(vi,V3, V4,V5 
(V2,V3, V4,V5 



(vi,V2, V3, V4,V5 



First, applying this procedure to the radical r we decompose the integration problem into two parts: the integration 
of the distribution with semisimple algebra g/r, then the integration of the restriction of distribution to the integral 
manifold with the solvable symmetry algebra r. 



The last step can be performed by quadratures. Moreover, every semisimple Lie algebra g/r is a direct sum of simple 
ones which are ideal in g/r. Thus, the Lie-Bianchi theorem reduces the integration problem to involutive distributions 
equipped with simple algebras of symmetries. Thus, integrating of system (2), become so much easy. 

The Lie algebra g is solvable and non-semisimple. It is solvable because if g*^^^ ~ (v^, [vi,Vj]) = [g,g], we have 
fl'"^'' = [flifl] = (vi, • • • , V5), and g(^^ = [g''"^\g'-"^-'] = (vi, V2,2v3), so, we have a chain of ideals g(^) D g'^^ D {0}, and 
it is non-semisimple because its killing form 

/O \ 





5 -i 
VO -8 17/ 



is degenerate. 

According to the table of the commutators, g has two abelian 2 and 3-dimensional subalgebras spanned by 
(vi,V2,V3) and (v4,V5) respectively such that the first one is an ideal in g, so, we can decompose g in the to 
semidirect sum of g = K'^ k M^ . 

7 Conclusion 

In this article group classification of Newtonian incompressible fluid's equations flow in turbulent boundary layers and 
the algebraic structure of the symmetry group is considered. Classification of r-dimensional subalgebra is determined 
by constructing r-dimensional optimal system. Some invariant objects are fined and the Lie algebra structure of 
symmetries is found. 



10 



References 

[I] M.L. Gandarias and M.S. Bruzon, Type II diddcn symmetries through weak symmetries for some wave equation, Communications 
in Nnonlinear Science and Numerical Simulation, 2009, article in press. 

[2] M. Nad,iafikhah and S.R. He,iazi, Symmetry analysis of cylindrical Laplace equation, Balkan journal of geometry and applications, 
2009. 

[3] M. Nadjafikhah, R. Bakhshandeh-Chamazkotia, and A. Mahdipour-Shirayeha, A symmetry classification for a class of (2+1)- 
nonlinear wave equation. Article in press, doi:10.1016/j.na.2009.03.087. 

[4] P.J. Olver, Equivalence, Invariant and Symmetry, Cambridge University Press, Cambridge University Press, Cambridge 1995. 

[5] P.J. Olver, Applications of Lie Groups to Differential equations. Second Edition, GTM, Vol. 107, Springer Verlage, New York, 1993. 

[6] L.V. OvsiANNlKOV, Group Analysis of Differential Equations, Academic Press, New York, 1982. 

[7] G.B. Whitham, Linear and Nonlinear Waves. New York: Wiley, p. 617, 1974. 

[8] D. Zwillinger, Handbook of Differential Equations, 3rd ed. Boston, MA: Academic Press, p. 132, 1997. 

[9] R. CiMPOlASU and R. Constantinescu, Lie symmetries and invariants for a 2D nonlinear heat equation. Nonlinear analysis, 68 
(2008) 2261-1168. 

[10] Q. Changzheng and Q. Huang, Symmetry reduction and exact solution of the affinc heat equation, J. Math. Appl. 346 (2008) 
521-530. 

[II] P.G.L. Leach, Symmetry and singularity properties of the generalised Kummer-Schwartz and related equations, J. Math. Appl. 348 
(2008) 467-493. 

[12] H. AzAD, M.T. Mustafa, Symmetry analysis of wave equation on sphere, J. Math. Appl. 333 (2007) 1180-1188. 

[13] E. Demetriou, N.M. Ivanova, C. Sophocleous, Group analysis of (2+1)- and (3+l)-dimensional diffusion-convection equations, 
J. Math. Appl. 348 (2008) 55-65. 

[14] A. KusHNER., V. Lychagin and v. Rubstov, Contact geometry and nonlinear differential equations, Cambridge University Press, 
Cambridge 2007. 

[15] www.wikipedia.com. 



11 



